Solitons in thermal media with periodic modulation of linear refractive index 
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We address the existence and properties of solitons in thermal media with periodic modulation 
of linear refractive index. Many kinds of solitons in such optical lattices, including symmetric and 
antisymmetric lattices, are found under different conditions. We study the influence of the refractive 
index difference between two different layers on solitons. It is also found that there do not exist 
cutoff value of propagation constant and soliton power for shifted lattice solitons. In addition, the 
solitons launched away from their stationary position may propagate without oscillation when the 
confinement from lattices is strong. 

PACS numbers: 42.65.Tg, 42.65.Jx 



I. INTRODUCTION 

Nonlocal solitons have been studied in many physical 
systems, such as photorefractive crystals [l|, nematic liq- 
uid crystals lead glasses 04^, Kerr- like media p^ . 
materials with quadratic nonlinearities |lll4l4ll, atomic 
vapors [l^, and Bose- Einstein condensates [la, 113 j etc. 
They show some special and novel properties both in 
experiment and theory, for instant, the large phase 
shift [13 , self-induced fractional Fourier transform [Toj . 
attraction between two dark solitons [13, |2l|, etc. Re- 
cently, nonlocal surface solitons propagating along the in- 
terface formed by a nonlocal nonlinear medium and a lin- 
ear medium are found and observed [23 - [25| . This kind of 
solitons presents strong asymmetry including their pro- 
files and nonlinear refractive indices because of the dif- 
ferent boundary conditions. They are stable only when 



24| 



the number of intensity peaks is less than three 

Solitons are also found in some media with periodic 
nonlinearity (26l - l30| . Recently, solitons in nonlocal non- 
linear lattices are studied, such as odd, even, dipole, 
tripole and twisted solitons (sTI - fs^ l. For lattice solitons 
formed in nonlocal media with an imprinted modulation 
of linear refractive index, there exists a cutoff point of 
propagation constant due to the existence of the band- 
gap structure of periodic lattice [3l| . In thermal lattices 
with periodic nonlinearity, there can be found asymmet- 
ric solitons, which is shifted solitons. There exist a cutoff 
value of propagation constant and a minimal power for 
this kind of solitons in such lattices (s^. It is found that, 
in thermal nonlinear lattices, solitons can exist anywhere 
inside a band and two families of solutions bifurcate from 
each band edge because of the infinite range of nonlocal- 
ity and the presence of boundaries [s^. At the same time, 
surface lattice solitons are found at the interface between 
a uniform medium and a nonlinear lattice or at the edge 
of layered nonlocal media [13, HI] . Such solitons present 
strongly asymmetric shapes and they are stable in large 
parts of their existence region. There is no restriction on 
the number of peaks in stable surface solitons localized 



at a thermally insulating interface between layered ther- 
mal media and a linear dielectric, which is quite different 
from that in uniform nonlocal media [24!|. 

Interface solitons propagating at the interface of one- 
dimensional thermal nonlinear media with a step in linear 
refractive index have been studied in [SG*! . It is note that 
interface solitons in such media are asymmetric because 
of the difference of the linear refractive index. For tripole 
and quadrupole interface solitons, there exist two differ- 
ent types of solutions under some special conditions . 

In this paper, we study the propagation of solitons in 
the semi-infinite band gap of thermal media with periodic 
distribution of linear refractive index. In our model, the 
period modulation function of linear index is a step func- 
tion. We find the odd, even, and shifted solitons in this 
kind of media. In addition, there do not exist cutoff value 
of propagation constant for shifted solitons, which is dif- 
ferent from that in thermal media with periodic nonlin- 
earity. We also study the influence of the index-difference 
of the lattices and the boundaries on solitons. 



II. THEORETICAL MODEL 

We consider a (l-l-l)-dimensional layered thermal sam- 
ple occupying the region —L < x < L imprinted a lin- 
ear refractive index lattice. A laser beam propagating 
along the z axis is governed by the nonlocal nonlinear 
Schrodinger equation (NNLSE) for the dimensionless am- 
plitude q of the light field, coupled to the equation for 
normalized nonlinear refractive index variation n. 
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where x and z stand for the transverse and longitudinal 
coordinates scaled to the beam width and the diffrac- 
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tion length, Ud represents the difference between the 
hnear refractive indices, S{x) is a periodic function for 
the modulation of the linear refractive index. In our 
model, we study the cases of S{x) — H[cos{ttx / d)] and 
d{x) = [sin(7rx/(i)], which present the symmetric and 
antisymmetric lattices, respectively. Here H[£] is the 
Heaviside step function defined as 



Hm = 




(3) 



d is the layer width and the lattice period is 2d. For 
simplicity and without losing of generality, we use inte- 
ger value of d in this paper, which possibly leads to the 
incomplete layers at the two boundaries. Anyway, the 
forth term of Eq. ([!]) represents the periodic modula- 
tion of two different linear refractive index in the sample. 
This kind of modulation is different from that in Ref. [33| . 
When 6{x) = 1, it represents the higher linear refractive 
index, while 6{x) = 0, it represents the lower linear re- 
fractive index. If Ud — 0, Eq. ([T|) represents the NNLSE 
in uniform thermal media. 

We assume the two boundaries are thermally conduc- 
tive and the laser beam is far away from the bound- 
aries, so the boundary conditions are q{±L) — and 
n{±L) — 0. We search for soliton solutions of Eq. ([T]) 
numerically in the form q{x,z) — u'(x) exp(i6z), where 
w{x) is a real function, h is the propagation constant. 
The iterative method is used to get numerical solutions 
for different n^, d and h as shown in Figs. [TJH To elu- 
cidate the linear stability, we search for the perturbed 
solution of Eq. ([T]) in the form q— (w + u-\-iv) exp{ibz), 
the real part u{x, z) and the imaginary part v(x, z) of the 
perturbation can grow with a complex rate a upon prop- 
agation. Substituting the perturbed solution into Eq. 
(H)), one can get the linear eigenvalue problem around 
stationary solution 'w{x). 



au = —773—2 + ~ nv ~ 5{x)ndV, 



1 d'^v 
2dx' 
1 d^u 
2d^ 



bu + nu + 6{x)ndU + wAn, 



) (4) 



where An = —2 J_j^G{x, x')w{x')u{x')dx' is the refrac- 
tive index perturbation, the response function G{x, x') — 
{x + L){x' - L)/{2L) for x < x' and G{x,x') = {x' + 
L){x — L)/{2L) for x>x'. Ur (real part of a) represents 
the instability growth rate. 



III. LATTICE SOLITONS 

First, we study the solitons in the case of b(x) = 
7f[cos(7rx/(i)] (symmetric lattices) and d = 2 as shown 
in Fig. [TJ When > the center layer is with higher 
refractive index, and the number of higher-index layers 
is odd. In Figs. [IJa) and[T](c), there exist odd solitons in 
the center of the sample because intensity peaks always 




FIG. 1: (Color online) Profiles of a fundamental soliton (a) 
at nti = 1, 6 = 0.5, even solitons (b) at Ud ~ —0.5, a tripole 
soliton (c) at ritj = 2, b = 1, and dipole solitons (d) at = 
—0.5. (e) Energy flows versus propagation constant of (a) 
(solid red line) and (c) (dashed blue line), (f) Real part of 
the perturbation growth rate of (b) (solid red line) and (d) 
(dashed blue line). In all cases d = 2. 



reside in the higher-index of the layers. Their profiles are 
similar to those in other types of optical lattices, while 
the distribution of nonlinear refractive indices are differ- 
ent [31I, H^l ■ When Ud < the center layer is with lower 
refractive index, and the number of higher-index layers is 
even. In Figs. \Vlh) and[ljd), one can find even solitons, 
including symmetric ones [Fig. [TJb)] and antisymmetric 
ones [Fig. [Ud)], at the center of the sample. As we know 
that the beam width decreases as the propagation con- 
stant increases. In optical lattices, the smaller the b is, 
the more layers the soliton profile occupies. That is the 
number of the intensity peaks decreases as the propaga- 
tion constant increases [Figs, mb) and[IJd)]. 

The linear stability analysis shows that fundamental 
and tripole lattice solitons in Figs. [Ua) and [IJc) are 
stable in their existence domain for > 0. Energy flows 
U — [(/l^dx of fundamental and tripole solitons are 
monotonically growing functions of b as shown in Fig. 
[TJe). However, when Ud < 0, the even and dipole solitons, 
as shown in Figs. [IJb) and[TJd), are unstable for Ud — 
—0.5 [Fig. [Ijf)]. We find that the antisymmetric dipole 
solitons as shown in Fig.[IJd) can be stable for Ud < —0.6 
due to the confinement of the deep lattice, whereas the 
even solitons as shown in Fig. [Hb) are unstable for any 
value of Ud- 

Then, we discuss the solitons formed in the antisym- 
metric lattices. The modulation of linear refractive index 
is governed by i5(a;) = H[sm{Trx/d)] and lattice structure 
is shown in Fig. [2] with d = 2. If there is only one pe- 
riod in the sample, i.e. d — L, the antisymmetric lattice 



FIG. 2; (Color online) Solitons in antisymmetric lattices at 
(a) Ud ~ —0.5, (b) Ud = —1, (c) and (d) b = 0.5. In all cases 
rf = 2 and i = 30. 



FIG. 3: (Color online) Propagations of solitons launched away 
from the sample center for (a) and (b) b = 0.5, = 1, d = 1, 
(c) b = 1, nd = 1, d — 2, and (d) 6 = 0.5, = 3, d = 3. In 
all cases L — 30. 



will be the same as the model in [36|. Here solitons are 
always asymmetric because of the antisymmetry of lat- 
tices. Figures [^a) and ^h) show the fundamental and 
even asymmetric solitons, respectively. One can see that 
the main intensity peaks reside in the higher-index lay- 
ers, and the larger the propagation constant the higher 
the degree of symmetry. Figures [5Jc) and [H^d) show the 
high-order asymmetric solitons at the same propagation 
constant. One can see that the degree of the restriction 
of the lattices is proportional to | | , and the asymmetry 
of the solitons is almost independent of rid- The funda- 
mental and the high-order asymmetric solitons shown in 
Figs, d^a), djc) and^d) are stable in their whole exis- 
tence domain, while the even solitons shown in Fig. [2{b) 
are unstable. 

As we know in [38], a soliton launched away from the 
sample center will oscillate periodically because of the 
restoring boundary force. In our model, we simulate 
the beam propagations with finite difference method in 
the symmetric lattices. Solitons launched away from the 
sample center undergo oscillation, as shown in Figs, ^a.) 
and [3]Jb) , when the boundary force is stronger than the 
constraint force of the higher-index waveguide 39] . It is 
important to note that the intensity peaks jump between 
the adjoining higher-index layers during propagations be- 
cause the peaks will not reside in the lower-index layers, 
which is the same as discrete diffraction i40j]. If the soli- 
ton is launched further away from the sample center as 
shown in Fig. [DJb), it also oscillates periodically with 
larger span, but its oscillating period is the same as that 
in Fig. ini^a). 

If we increase the layer width [Fig. [Sjc)], or the differ- 
ence between the two linear refractive indices [Fig. [31[d)], 
the solitons willpropagate without oscillation. Similar to 
the results in [3^, when the constraint force is stronger 
than the restoring force, the intensity peak is difficult to 
jump into other higher-index layers for large and d. 
It is worth to note that only fundamental solitons will 
oscillate periodically, the high-order solitons will oscil- 



late irregularly and destroy themselves for small Ud and 
d. However, stationary high-order solitons can form for 
large Ud and d when they are launched away from the 
sample center. 

Based on the above analysis, when the constraint force 
of the lattices is large, solitons can form asymmetri- 
cally away from the sample center, that is shifted soli- 
tons. Figures lUJa) and lU^b) show the profiles of funda- 
mental, dipole, tripole, and quadrupole shifted solitons 
at the same conditions. Their mass center, defined as 
Xm = /.!^2:|gpda;/ \q\'^dx, are 3.1, 6.8, 3.9, and 7, 
respectively. This kind of solitons exists when the con- 
straint force of the higher-index waveguide is stronger 
than the restoring force exerted by two boundaries due 
to the shift of the soliton center. Figure Sljc) shows that 
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FIG. 4: (Color online) Profiles of (a) fundamental and tripole 
with b — 0.05, Ud = 2, and d = 3; (b) dipole and quadrupole 
shifted solitons with 6 = 0.05, Ud = —2, and d = 3. (c) 
Energy flows versus propagation constant of fundamental, 
dipole, tripole, and quadrupole solitons. (d) Real part of the 
perturbation growth rate of the quadrupole soliton. 
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FIG. 5: (Color online) Profiles of shifted solitons at (a) b = 0.5 (black line) and 6 = 1 (red line), (b) b = 0.5 (black line) and 
6=1 (red line), and (c) b = 0.5. In all cases = 5. (d)-(f) Propagations of shifted solitons corresponding to (a)-(c) at b = 0.5. 



their energy flows are monotonically growing function of 
b and there do not exist a cutoff value of b for any shifted 
solitons. It is found that there exist a minimal power 
and a cutoff value of b for shifted solitons in the model 
in (3^ , because the modulation in their model is induced 
by the nonlinearity, which will become weak for small 
propagation constant and low soliton power. However, 
in our model, the modulation of linear refractive index 
does not depend on the constant propagation, so there 
does not exist a cutoff value of b for shifted solitons. For 
fundamental, dipole, and tripole shifted solitons shown 
in Figs, mja) and|Ub), they are all stable in their whole 
existence domain, while for the quadrupole shifted soli- 
ton, there exists an instability region shown in Fig. Hl^d). 
We can see that the instability region is dependent on 
rid, and the instability region decreases as increasing of 
\nd\- As |nd| > 8, almost all solitons are stable. 

Another type of shifted solitons formed away from the 
sample center is shown in Fig. [5] They are separated by 
the central layers and located in any higher-index layer, 
including symmetrical [the red lines in Figs. [SJa) and 
[5]Jb)], antisymmetrical [the black lines in Figs. [Slja)and 
[5l^b)], and asymmetrical [Fig. [^c)] ones, because of the 
restraint force of lattices and the existence of boundaries. 
Almost all intensity peaks reside in the higher-index lay- 
ers for this type of solitons, and they do not affect other 
peaks during the propagation due to the restriction of the 
lattices. For this reason, the solitons shown in Figs. EJa)- 
[5]Jc) are stable including the quadrupole shifted solitons, 
which is different form that in Fig. IDJb). Their propaga- 
tions are shown in Figs. [5][d)-[5][f). One can see from Fig. 
[5]Jf) that the two intensity peaks are stable and have no 
impact on each other, even if they are far away from the 
sample center. 

Finally, we discuss the influence of the index-difference 
on solitons with d and b fixed in the antisymmetrical lat- 
tices. When Ud is small, as shown in Fig. IHJa), there 
are three main intensity peaks with a significant part of 
the field residing in the lower-index layers and exist sev- 



eral sub-peaks. That is because the higher-index layers 
present a little restriction at small Ud- As nd increases 
[Fig. ini[b)], the part of the field residing in the lower- 
index layers decreases, and the middle intensity peak 
decreases gradually. When Ud = 10, almost all soliton 
power resides in the higher-index layer and the middle 
peak disappears as shown in Fig. ^c). It shows that 
the restriction of the higher-index layers increases as Ud 
increases. Figure HJd) shows the relationship between 
nd and beam width and the energy flow. One can see 
that the energy flow is proportional to Ud, while the 
beam width changes non-monotonically with increasing 
Ud- When Ud < 1.5, the beam width decreases monoton- 
ically as Ud increases as shown in Fig. [SJ^d). The main 
reason is that the widthes of two lateral peaks decrease 
as Ud increases, while the width of the middle peak is 
unchanged [Fig. [5l[a)]. When Ud > 1.5, the beam width 
increases monotonically as Ud increases [Fig. [SJ^d)] be- 
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FIG. 6: (Color online) Profiles of solitons at (a) Ud = 0.5 and 
Ud = 1.5, (b) Tid = 3 and Ud = 7, and (c) Ud = 10. In all 
cases d — 2 and b = 0.5. (d) Energy flow (blue dashed line) 
and beam width (red solid line) versus n^. 
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cause the width of the middle peak decreases gradually, 
while the widthes of two lateral peaks are unchanged with 
increasing Ud- When > 5, the beam width approaches 
5 because the middle peak almost disappears and there 
is little change of the soliton profile any more [Fig. [5{b) 
and Fig. El^c)]. For comparison, soliton profiles in sym- 
metric lattices are symmetric when the constraint force 
is weaker than the restoring force [Figs. [T], whereas they 
can be either symmetric [Fig. [5] or asymmetric [Figs, 
m and [5] when the constraint force is stronger than the 
restoring boundary force for large Ud [Fig. |4]. For soli- 
tons in antisymmetric lattices, however, they are always 
asymmetric [Figs. [2] and [6]. 

IV. CONCLUSION 

To conclude, we have studied the propagation of soli- 
tons in thermal media with periodic distribution of linear 
refractive index including symmetric and antisymmetric 



lattices. In this kind of optical lattice, we found fun- 
damental, odd, even, and shifted solitons. The shifted 
solitons are found do not exist cutoff value of propaga- 
tion constant and power, which is different from that 
in thermal lattices with periodic nonlinearity. We also 
studied the solitons launched away from their stationary 
position. Their propagation is related to the difference 
of the refractive indices. 
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